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Q"^ ■ A explicit formula on semiclassical Green functions in mixed position and 

momentum spaces is given based on Maslov's multi-dimensional semiclassi- 
0-i. cal theory. The general formula includes both coordinate and momentum 
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Semiclassical methods are as old as quantum mechanics and an enormous 
amount of work has been done in the field. Among them, the WKB method 
can be found in any standard quantum mechauics textboob.B But it IS usu- 
ally confined in one dimensional case. There has long been a need for the de- 
velopment of a method that can be used to deal with multi-dimensional prob- 
lems. In last two decades, the Maslov and Fedoriuk's (MF) n-dimensional 
semiclassical theorvn has achieved a remarkable success in atomic and molec- 
ular physics field.H The semiclassical method emphasizes a suitable choice 
of representation space (usually a mixed position-momentum space) to over- 
come the semiclassical diversity in configuration or momentum space. Delos 
has shown how to construct a uniform semiclassical wave function by MF's 
approach.^ The semiclassical Green functions or propagators in configura- 
tion and momentum space has been available for a long time.! But they 
can not work in the singular regions such as near caustics or foci of classical 
trajectories. Semiclassical Green functions and trace formulas in connection 
with caustics are analyzed in detail by Littlejohn and coworkersi.0 They em- 
phasize geometrical properties of the semiclassical trace formula and shows 
the intrinsic property of Maslov index. In this paper we want to construct a 
semiclassical Green function explicitly in mixed position-momentum spaces, 
which can be used locally to fix the diversity of its counterpart in config- 
uration or momentum space at any boundaries between classically allowed 
and classically forbidden regions. The Green function in mixed space is im- 
portant to construct to a uniform semiclassical Green function in customary 
configuration space. In fact, in above mentioned singular regions we can 
obtain a uniform semiclassical Green function in following steps: Firstly, we 
transform the primitive Green function in configuration space into a suit- 
able chosen mixed space. Then a finite Green function in the mixed space 
is calculated. Finally, we transform it back into configuration space, and 
sewing it with the primitive Green function smoothly to construct a uniform 
semiclassical Green function. In the following we will achieve the three steps 
explicitly. 

For simplicity's sake we consider a single particle without spin in n- 
dimensional configuration space, and its Hamiltonian independent of time. 
The propagator at constant energy or Green function G {q', q", E) satisfies 
the inhomogeneous equation 

[E-Hip",q")]Giq',q",E) = 5iq"-q') (1) 

where q' is the initial position, q" the final position and E is the total energy. 
The semiclassical approximation G of G is given byi 
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Giq",q',E) 
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(n+l)/2 



class, traj. 



(2) 



where 



Siq",q',E) 



pdq, 



(3) 



is the classical action evaluated along the classical path which leads from q' 
to q" at the given energy H{p, q) = E and /i is the Maslov index, which is 
related with topological property of the ray of trajectories. Its geometrical 
properties has been discussed in detail by Littlejohn and coworkerji.0 The 
summation in Eq.(|^) is over all classical trajectories from q' to q" . The 
determinant is defined as 



dq'dq" 
dEdq" 



dq'dE 



(4) 



which is a n + 1 th determinant and can be considerably simplified as n — 1th 
one if introducing a local coordinate system centered on the trajectoryffl.i 
Denote the coordinates (y, z) , with z runs along the trajectory and n — 1 
coordinates y transverse to the path in such a way y = specifies the path. 
In the local coordinate, 
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z'z" 



dy'dy" ' 

In fact, above restrictions on the local coordinate are not necessary and have 
been shown in Refs.(6) and (7). 

The semiclassical propagator in momentum space is given , in analogy of 
its counterpart in configuration space, by 



F{p",p',E) 
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where 



T{p",p',E) 



qdp, 



(7) 



is the classical action along the classical path from p' to p" at the given 
energy H{p, q) = E and u is the Maslov index calculated in momentum 
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space. Usually, z/ 7^ /i because of the different topological properties of 
trajectories between configuration and momentum spaces. The determinant 
Dt is defined in an analogy way as 



dp' dp" 
dEdp" 



dp'dE 



Applying a similar local coordinate in momentum space, say, Pz and Py are 
along and transverse the trajectory respectively, we have the stability matrix 
in momentum space 



WzPzl 



(9) 



According quantum mechanics Green functions both in position and mo- 
mentum spaces must be finite. However, semiclassical approximation of them 
contain divergences at any boundary between classically allowed and classi- 
cally forbidden regions, and at any caustics or foci of classical trajectories. 
The reason causing such divergences is not the semiclassical theory itself, but 
configuration or momentum representations are not suitable in these singu- 
lar regions. If we transform the problem to some mixed position-momentum 
representation, such divergences can be repaired. It has been shown that 
there "almost always" exists a representation in which the semiclassical ap- 
proximation does not divergent in these singular regionJi.i 

We choose mixed coordinates {pa, qp), where a = 1, - ■ ■ ,k and (3 = k + 
1, • • • ,n,i.e, which are disjoints set. So the set PaQ/s contains no canonical 
conjugate pairs. The propagator or Green function in the mixed space can 
be acquired from the counterpart in configuration space by partial Fourier 
transformation. The propagator F [p'l^q'p^p'^q'pE^ from p'^q^ to p'^q^ while 
propagating with the energy is defined by 



F{p%M,E) = {27^%)-' j dq^j dq'^G{q"q'E) (10) 

xexp [i {p'Jo.-p'Lql) /h]. 

Its semiclassical approximation F {^'^q'p^p'^q'pE^ can be acquired by inserting 

G {q"q'E) into Eq.(|l^) and evaluating the integral by A;— dimensional station- 
ary phase approximation(SPA). Firstly, we calculate the integral on q'^ , that 
is 



F {q'\pWpE) ^ i2nh)-'/' J dq'^G iq"q'E) exp {p'J^ - p'^) ^ • (11) 
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where F {^q" .p'^q'^E^ represents the integral after SPA. Inserting Eq.(0) into 
the expression, the exponent becomes 



~S = S{qyE)+p'J^-p'^ql 



(12) 



apart from the factor i/h. It is noted that we adopt a single S to denote 
the exponent before and after SPA. For convenience, Maslov index and the 
summation over classical orbits have been dropped temporarily. The SP 
points g^Q satisfy 



dS 



dq'a 



-Pa 



a 



k. 



After SPA, the preexponential factor can be written as 
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{27rh) 
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A/2 



1/2 



{2mhf^-''^/''"'' V dp',, 
To simplify the expression further, we rewrite as 
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and the second derivatives of S now should be replaced the derivatives of S. 
To achieve the replacement, we can use the following relations 
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, Eq.(pl) 


becomes 
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dq'^dq'; 
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dq'^dE 
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iD'.f" (17) 



It is noted that in the first fine S is replaced by S . After SPA, the integral 
in Eq.(pJ]) becomes 

where 



S = S{q"q^QqpE) + pj^^ - fj'^. (19) 

is the value of Eq.(|l^) at the first SP point. Secondly, we again approximate 
the integral over by SPA to get the mixed Green function 



(20) 

Again F [p'^q'^^p'^q'^E^ denotes the integral after SPA. The SP points now 
satisfy 



dS 



dq'L 



-P'L: a = l,---,k. 
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The relations in Eq.(|TB|) are replaced by 
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(22) 



which leads to the replacement of the derivatives with respect to q'^ in the 
first column of D'^ by the ones to p'^. Then the preexponential factor arrives 
at 
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It is noted that the factor i'^ after two integrations has been absorbed in the 
final determinant . Therefore, the final Green function in PaQ/i space is 
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where 



S = S {q'^o<lW»o<l'i3E) + pWao - P'Wao- (25) 
is the value of Eq.(|r^) at the second SP point though with the same symbol. 



In Eq.(p^) the Maslov index u has been added and is related with the topo- 
logical properties of classical trajectories in the mixed space. The summation 
over classical trajectories has also been attached. For convenience we have 
adopted the single symbol S to represent the phase factor before and after 
SPA. S can be written in a regular form 



S 



9/3 



Qadpa 



(26) 



which is the classical action evaluated the trajectory from p'^q'/^ to p^g^ at 
given energy H {pq) = E. It can be acquired by differentiate Eq.(^) with 
respect to the terminal coordinates. Eq.(P3D can be rewritten as 



S 



ni/s 
Padqa+ / Ppdq/s + p'^q'^Q - p'X' 



(27) 



which has contained the first term in Eg . (|2q) . Differentiating it with respect 
to p'^ , we have 



dp'^ 

Similarly, one finds that 



Pet a„/ ^ 'iaO ^ Pa 



dq'a 
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dp' 
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(29) 
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Eq.(0) and (|2|) immediately leads to (|2|). 

The stability matrix Dg can also be simplified Ith determinant 

depending on the choice of local coordinate centered on the trajectory in the 
mixed space. If we choose a position coordinate z along with the trajectory 
in the mixed space, {pa-, y/3) transverse the path, then 



Ds 



1 



\z'z"\ 



92 S 



925 



(30) 

Conversely, choosing a momentum coordinate along with the trajectory 



and {pya,(l/3) transverse the path, we have 



\P'zPz\ 



925 



dPyadp'yc, 
925 



dPyadq'^ 
925 



Formula (p4 ) is a unified expression. If a = 0, /5 = 1 ■ ■ ■ n, it will recover as 
Eq.(^, the Green function in configuration space. Conversely, if a = 1 ■ • - n, 
/3 = 0, Eq.(^), the propagator in momentum space will immediately emerge. 

For most semiclassical physical problems, it is convenient to work in con- 
figuration space. Eq.@ usually works quite will in most regions. But in 
some regions, such as near caustics and foci of classical trajectories, it will be 
divergent. In these regions, we can " almost always" fflchoose a suitable set 
of mixed position-momentum coordinates to calculate the mixed Green func- 
tion, which is not divergent. From the function, an accurate configuration- 
space Green function can be constructed be inverse partial Fourier transfor- 
mation, 



G{4WE) = {2nh)-'Jdp:Jdp'J{pW;,,p'J,E) (31) 

Xexp {-pWa + P'Wa) M ■ 

Then we can sew it with the primitive Green function smoothly to get a 
uniform Green function. 

The author thanks Professor M. L. Du for helpful discussions. 
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